Linking high and low temperature plasticity in bulk metallic glasses: thermal 
activation, extreme value statistics and kinetic freezing 
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At temperatures well below their glass transition, the deformation properties of bulk metallic 
glasses are characterised by a sharp transition from elasticity to plasticity, a reproducible yield 
stress, and an approximately linear decrease of this stress with increasing temperature. In the 
present work it shown that when the well known properties of the under-cooled liquid regime, in 
terms of the underlying potential energy landscape, are assumed to be also valid at low temperature, 
a simple thermal activation model is able to reproduce the observed onset of macro-scopic yield. At 
these temperatures, the thermal accessibility of the complex potential energy landscape is drasti- 
cally reduced, and the statistics of extreme value and the phenomenon of kinetic freezing become 
important, affecting the spatial heterogeneity of the irreversible structural transitions mediating 
the elastic-to-plastic transition. As the temperature increases and approaches the glass transition 
temperature, the theory is able to smoothly transit to the high temperature deformation regime 
where plasticity is known to be well described by thermally activated viscoplastic models. 

PACS numbers: 62.20.D,81.05.Kf, 46.25.Cc 



I. INTRODUCTION 

The low temperature mechanical properties of struc- 
tural glasses, such as Bulk Metallic Glasses (BMGs), are 
characterised by a large elastic strain and a well defined 
yield point^. This latter feature is particularly the case 
for tensile deformation, but work has also shown that 
if a pure uni-axial stress state can be achieved in com- 
pression, the yield point is also reproducible with a high 
WeibuU modulus comparable to that of a work-hardened 
crystalline metaP. Thus, although BMGs lack good duc- 
tility, they are not brittle in the ceramic sense since the 
relatively low WeibuU modulus of a ceramic indicates a 
failure mechanism sensitive to macroscopic flaws. Like 
work hardened crystalline metals, this suggests that the 
structure of a BMG is relatively homogeneous above a 
certain length scale, and yield is an intrinsic material 
property. Other related properties characterising the low 
temperature deformation regime of BMGs are a strong 
spatially heterogeneous plasticity, a weak strain rate sen- 
sitivity of the yield stress and its approximately linear 
decrease with respect to increasing temperature. 

There has been much work attempting to characterise 
the underlying microscopic structural mechanisms medi- 
ating plastic evolution in BMGs. In the seminal work 
of Spaepen^, the picture of single atom migration into 
regions of excess free volume was proposed. In subse- 
quent theoretical work by Argon^, the atomic scale struc- 
tural transformation was hypothesised to consist of a 
rearrangement of atoms which mainly resulted in a lo- 
cal shear stress relief. The early atomistic work of Talk 
and Langer— supported this picture introducing the con- 
cept and terminology of the Shear Transformation Zone 



(STZ), in which groups of atoms collectively undergo a 
local shear transformation. Static and ultra-high strain 
rate dynamic atomistic simulations have given much in- 
sight into the concept of the STZ — how it may be re- 
alised at the atomic level and how together they might 
collectively behave for the emergence of material yield 
and subsequent plastic flow^"— . Indeed, recent molecu- 
lar dynamics simulation work by Guan et alJ^ suggests 
that at high enough strain rates macroscopic yield cor- 
responds to a stress-induced reduction of the glass tran- 
sition temperature, providing a link to the quite general 
phenomenon of jamming/unjammingi^^. 

In both the works of SpaepenS. and Argon^, the theory 
of thermal activation played a central role, an approach 
that forms the basis of many numerical mesoscopic mod- 
els of quasi-static plastic flow in BMGs^"— . Indeed for 
temperatures outside the low temperature deformation 
regime, and close to the glass transition, simple ther- 
mally activated viscoplastic models can describe well the 
experimentally observed strong temperature and strain 
rate dependent deformation behaviouri^iiS^ — a regime 
of deformation that is now spatially homogeneous and 
closely related to the viscosity of the under-cooled liquid 
regime. 

In the under-cooled liquid regime, two distinct relax- 
ation timescales are believed to determine the physics 
of viscosity. These are the slow a-relaxation processes 
and the fast /3-relaxation processes, both of which were 
first postulated by considering the underlying potential 
energy landscape (PEL)^'^, and latter confirmed by ex- 

The PEL viewpoint has afforded great physical insight 
into the origin of the strong temperature dependence of 
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viscosity in approaching the glass transition temperature 
from the under-cooled liquid regime^iiS^, as evidenced 
by the developed understanding of the well known An- 
gell plot that introduced the important material param- 
eter of fragility^^ . The PEL for an under-cooled liquid is 
seen as a complex energy landscape of hills and valleys 
consisting of mega-basins, whose exit characterises the 
a-relaxation, and whose underlying finer structure char- 
acterises the /3-relaxation2^i2ii2^i2£. It is in this context 
that /3-relaxation is said to mediate a-relaxation, a view- 
point that is also extended to temperatures below, but 
close to, the glass transition. In this high temperature 
regime (of a BMG), /3-relaxation is seen as being equiva- 
lent to STZ activity causing a build up of internal stress 
that ultimately underlies the a-relaxation energy land- 
scape^-. Indeed, for a wide range of BMGs, there exists 
a strong correlation between the apparent barrier energy 
corresponding to a-relaxation (derived from the exper- 
imental kinetic fragility) and the characteristic barrier 
energy associated with thermally activated plastic flow 
(derived from temperature dependent strain rate defor- 
mation behaviour 

As mentioned in the introductory paragraph, at tem- 
peratures well below the glass transition, distinctly dif- 
ferent deformation behaviour is observed. Because of the 
experimentally observed weak temperature and strain 
rate dependence of the yield (and maximum flow) stress, 
there has been some focus on athermal explanations^^i^^ 
in which an effective temperature is introduced arising 
from the inherent quenched structural disorder of the 
amorphous solid. However, there are experimental in- 
dications that a thermal activation mechanism also un- 
derlies the strongly inhomogeneous plasticity character- 
istic of this low temperature regime. This is particularly 
the case for compression geometries in which the initia- 
tion and propagation of shear bands and serrated versus 
non-serrated flow has been studied^"—. Thermal acti- 
vation theories have also been used to describe the weak 
reduction in yield stress with increasing temperature. Us- 
ing elasticity theory, Argon^ obtained a stress dependent 
critical barrier energy whose leading order term resulted 
in the reduction of yield stress following a T'^l'^ power- 
law^. More recently, Johnson and Samwer^^: have used a 
thermal activation model to derive an expression for the 
characteristic stress of plastic flow, which they argued is 
representative of the characteristic yield stress, to obtain 
a (T/Tg)^/'^ temperature dependence that describes well 
the universal trend seen in BMGs between yield stress 
(divided by a representative elastic modulus) and tem- 
perature. Again, the concept of a critical STZ barrier 
energy was used. 

In this paper, the properties of the under-cooled liquid 
PEL are further exploited to develop a theory of ther- 
mally activated plasticity for BMGs which is valid at both 
high and low temperatures below the glass transition. 
The current work extends on previous work^ and gives 
more detail to that presented in ref4^. Specifically, the 
complexity of the underlying PEL, in terms of the num- 



ber of underlying available irreversible structural trans- 
formations and their energy distribution, is extended to 
temperatures below the glass transition and used to pos- 
tulate a coarse grained distribution of a-relaxation bar- 
rier energies. When used in conjunction with the thermal 
activation hypothesis, a plastic transition rate is obtained 
that characterises a temperature and stress scale at which 
non-negligible plasticity occurs giving a measure of the 
yield stress and a definition to the critical barrier energy 
alluded to by Argon* and Johnson and Samwer— . Two 
distinct temperature regimes emerge, a high tempera- 
ture regime close to the glass transition temperature in 
which the yield stress rises rapidly and non-linearly from 
zero with decreasing temperature, and a low temperature 
regime in which the plastic transition stress increases ap- 
proximately linearly with decreasing temperature. It is 
found that this latter deformation regime is strongly af- 
fected by the universal phenomena of freezing and ex- 
treme value statistics of the PEL, leading naturally to 
an origin for the experimentally observed universal tem- 
perature dependent deformation behaviour of amorphous 
solids at low temperature. 



II. MODEL DEVELOPMENT 

A. Definitions and the development of a statistical 
framework 

To study the deformation properties of a structural 
glass, the characteristic time scale, Tp, associated with 
plastic activity at a particular length-scale is considered. 
For temperatures close to and above the glass transition 
temperature, the viscosity is proportional to such a time 
scale^^i^S. At temperatures below the glass transition 
temperature, where the initial plasticity is characterised 
by a transition from elasticity to a plastic flow regime, 
the inverse [rp] ^ can be viewed as a plastic rate. Such a 
basic description of low temperature plasticity cannot in- 
clude elastic/plastic accommodation issues and any cor- 
responding structural evolution, and is therefore unable 
to give quantitative insight into the nature of subsequent 
macroscopic plastic flow such as shear banding under 
compression and brittle fracture under tension. However, 
as a quantitative measure of the stress at which signif- 
icant plasticity begins to occur, such a description will 
suffice and indeed has been used in the past to describe 
the transient yield strength prior to macro-plasticity^^ 
and the yield rate by Johnson and Samwer—. 

If there exists a certain length scale above which a 
glass maybe considered homogeneous, then all material 
volume elements of a size comparable to this length scale 
are equivalent. It is the calculation of [tp^rve] ^ for such 
a representative volume element (RVE) that is of current 
interest. A RVE has well defined properties because a 
sufficient level of self averaging occurs with respect to a 
shorter length-scale which characterises the underlying 
heterogeneity. Let the volume of such a heterogeneity 
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length-scale contain a characteristic number of atoms, N. 
It is with respect to this heterogeneous volume element 
that the self averaging is performed, giving 

N' 

[rp,RVE]-' (T) = ^ [rp,„]-i (T) = N' x [r^]-' (T), (1) 

n=l 

where N' is a sufficiently large number of heterogeneous 
volume elements such that the average 

N' 

[^p]"'(7^) = ]^E[^p.«r'(T), (2) 

n=l 

has converged. Thus the RVE contains N' x N atoms. In 
eqns. [1] and [21 [Tp,n]~^ (T) represents the plastic rate of 
the nth heterogeneous volume element within the RVE. 
In the remainder of this paper it is [rp]"^ (T) that will 
be considered further. 

Under the assumption of thermally activated plasticity, 
the plastic rate associated with one particular realisation 
of a heterogeneity volume may be written as a linear sum 
of the M — M{N) individual irreversible transition rates 
available to that volume element: 

[-P.n]-^ (T) ^ f: [rpo.„.]-^ (T) exp (-^) , (3) 

where [Tpo,ni] ^ {T) and EpQ^ni are the attempt rate and 
barrier energy for the ith irreversible structural trans- 
formation within the nth heterogeneous volume element, 
respectively. Substitution of eqn. |3] into [2] then gives 

[rp]-i (T) = M ^[Tpo]-^ (T) exp (-^) ) , (4) 

where (• • • ) represents an average with respect to sam- 
pled realisations of heterogeneous volume elements. Such 
an average will be performed under the assumption that 
the heterogeneous volume elements are statistically inde- 
pendent from each other. 

To further simplify eqn. SI knowledge of the underlying 
PEL of an amorphous solid is required. Fig. [T] displays a 
typical one dimensional schematic of a microscopic PEL 
indicating the two characteristic energy landscape scales 
of the a-relaxation and /3-relaxation processes. This fig- 
ure schematically shows that to escape the mega-basins 
associated with the a-relaxation processes, multiple (and 
reversible) /3-relaxation (STZ) activty will occur^". Ex- 
iting such mega-basins results in irreversible structural 
transformations and thus emergent meso/macro-scopic 
plasticity. It is such structural transformations which 
enter explicitly into eqn. [3l and thus the barrier energies 
Epo^ni will represent those of the available a-relaxation 

processes (for example, the barrier energies 

and -^3""* schematically shown in fig. [T]), whereas the 
multiple mediating /3-relaxation/STZ activity will enter 
implicitly into the temperature dependent prefactors of 




reaction coordinate 



FIG. 1. Schematic of a one-dimensional realisation of the po- 
tential energy landscape of a structural glass. The horizontal 
reaction coordinate represents a trajectory in 3A^ dimensional 
space in which the system moves from minimum to minimum 
via a saddle-point configuration. The red-line represents the 
actual microscopic landscape and demonstrates that many 
intermediate minima (/3-relaxation/STZ processes) must be 
traversed before exiting a mega-basin (a-relaxation process). 
The blue dashed line indicates a coarse graining of the land- 
scape, in which only the a-relaxation landscape is retained. 



eqn. [3] This separation of energy scales and dynamics 
represents a natural coarse graining of the PEL in which 
the microscopic landscape of the /3-relaxation is effec- 
tively integrated out resulting in prefactors which are 
themselves diffusive and thus based on thermal activa- 
tion. 

The average in eqn. S] is with respect to the 
different possible temperature dependent pre-factors, 
[''po.ni] ^ [T) and mega-basin (a-relaxation) barrier en- 
ergies EpQ^ni and can be performed via a corresponding 
probability distribution with respect to these variables. 
Such a distribution naturally acknowledges that a BMG 
can admit a wide range of irreversible microscopic pro- 
cesses, a feature that has been discussed and exploited 
in past work on the early stages of plasticity^i^"— and 
more structural relaxation^"—. 

What form should this distribution take? Because the 
attempt rate for an identified escape path in the coarse 
grained landscape will depend on the specific geometry 
of the underlying microscopic landscape, its actual value 
and temperature dependence will be quite unique to that 
escape route and no one particular value is expected to 
be a representative case. This realisation results in the 
considerable simplification of allowing independent aver- 
aging of the pre-factor and the barrier energy, an ap- 
proximation that has also been rationalized by Gibbs 
et and whose consequences have been discussed by 
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Primak^*^. Eqn. |4] then becomes 



[.,P=M(M"^(r))(exp(-g,) 



(5) 



Since thermal activation of muhiple /3-relaxation/STZ 
processes mediates the a-relaxation process, a natm'al 
choice for the average prefactor will be 



[TpooJ exp 



pOO 



(6) 



where ii^poo is a characteristic barrier energy for j3- 
relaxation. This form is motivated by experiment which 
demonstrates that the average temperature dependence 
of /3 relaxation is Arrhenius both above and below the 
glass transition temperature2^i21i^. The lack of a tem- 
perature dependence in the corresponding characteristic 
pre-factor is also motivated by the same assertion lead- 
ing to the independent averaging within eqn. 31 where no 
particular pre-factor (and its temperature dependence) 
is representative — a situation that upon averaging nat- 
urally leads to little or no systematic temperature de- 
pendence. It will turn out, that any temperature depen- 
dence that might arise from the averaging of eqn. |6] will 
play only a minor role in the underlying properties of the 
model. 

Thus the final equation for the average plastic rate of 
a heterogeneous volume is 



[TpooJ exp 



-5'poo \ 
kBTj 



M ^exp 



_Epo \ 
kBT) 



^ (7) 

consisting of a diffusive attempt rate with a simple Ar- 
rhenius temperature dependence, representing the medi- 
ating /3-relaxation dynamics, and a thermal factor whose 
temperature dependence will be derived from the statis- 
tical properties of the a-relaxation coarse grained PEL. 



B. A Gaussian distribution of barrier energies, 
extreme value statistics and the emergence of 
kinetic freezing 

How is M related to TV and what is the form of P{E)1 
There will exist a relation between M and N whose 
origin lies in the properties of the PEL. Stillinger has 
shown the number of minima (inherent structures) that 
a structural glass can admit scales exponentially with 
system size^""— . Later work by Scott Shell et ex- 
tends such exponential scaling to saddle points, as was 
also done for a many-particle random Gaussian potential 
by FyodoroA*^. A toy model that approximates the min- 
ima of the PEL as an N hyper-dimensional cube^ has 
also been extended to the enumeration of saddle points 
resulting in a similar exponential scaling of the number 
of saddle points'^^'. Thus the relation between M and N 
is generally excepted to be exponential, M = exp{aN), 
where for the case of counting minima (inherent struc- 
tures), a is a measure of the configurational entropy per 



atom. Such an exponential scaling of the number of avail- 
able transition paths is not new and has been assumed in 
thermodynamical treatments of the viscosity of an under- 
cooled liquid^^'^®. In the present case, M counts the 
number of first order saddle points in the coarse grained 
PEL and a can be referred to as the configurational bar- 
rier entropy per atom. 

The form of P{E) turns out to be intimately related 
to the above exponential scaling, where the number den- 
sity with respect to barrier energy of a-relaxation barrier 
energies may be written as M{E) — MP{E). Similar to 
Stillinger 's work on the distribution of inherent structures 
as a function of system energy^, the work of Scott Shell 
et a^, and Derlet and Maass^ (within the frame work of 
a toy model) , have shown that such a number density can 
be well approximated by a Gaussian distribution, giving 



PiE) 



1 



: exp 



{E-{E)f 
2{E^} 



(8) 



with its first and second cumulants being extensive quan- 
tities, giving the mean as EN and the standard devi- 
ation as SE^/N. Here E can be viewed as the mean 
a-relaxation barrier energy per atom. Whilst such a dis- 
tribution will describe well the true distribution in the 
region where it is non negligible in value, it will give an 
overestimate of the barrier probability in the low barrier 
regime, indeed at a barrier energy equal to zero it will 
give a non-zero probability. Whilst this aspect will be 
addressed in sec III CI further development will assume 
that EN and SEVn have values such that P{E < 0) 
is negligible. In any event, the defining property of the 
distribution is not its precise form, but rather that it has 
an extensive first and second cumulant. 

Inspection of eqn. [7] reveals that the plastic rate is pro- 
portional to A/(exp(— /?£')) (here /3 = l/{kBT)) where 
(cxp(— /?£')) is the generating function of the distribu- 
tion evaluated at — /?. For the Gaussian distribution the 
generating function has a closed form, resulting in 




= M 



dEP{E)cxp(^- 



E 
kWT 



E 
kWT 



5E^' 



1 

2 \kBT 



(9) 



In the above the approximation arises due to the lower 
integral limit being extended from zero to — oo. 

For a sufficiently large N , the above (and [t^] ^ via 
eqnl?]) will have a negligible value when the factor within 
the argument of the exponential is negative and a large 
value when it is positive. Assuming E, SE and a are fixed 
material parameters, the critical temperature at which 
this occurs is defined by the condition 



E 



1 f SE 



fceTc 2 \kBTc 



-0, 



(10) 
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FIG. 2. Plot of a) plastic and freezing transition temperature, 
and b) apparent barrier energy evaluated at both the plastic 
and freezing transition temperatures as a function of a. 



whose relevant solution is, 



E 



{5E/'E) 



l-Jl-2 (SE/E) a 



(11) 



The critical temperature is thus independent of N and 
therefore the size of the heterogeneous volume element. 
Tc will be referred to as the plastic transition temper- 
ature. Fig. [2^ plots eqn. [11] as a function of a for the 
special case oi E — SE. For a — the plastic transition 
temperature diverges and at its maximum allowed value 
of one half, Tc — E/k^. Thus with an increasing a (an in- 
creasing number of available structural transformations) , 
Tc reduces. 

This latter result may be understood by investigat- 
ing the first and second moment of the distribution, 
P{E) exp(— which is straightforwardly given by the 
general expression for the nth moment. 



E 



(") 



1 



5"G(-/3) 



(12) 



where G{—/3) — (exp(— /3£^))) is the generating function 
for the Gaussian distribution and /3 — l/(fcBT). This 
gives the first moment as 



E^'\T) = N (e - 



and the second moment as 

E(^\T) = NSe'^ + Ie 



SE 



(13) 



(14) 



The first moment will be referred to as the apparent bar- 



rier energy. 



S(i)(T) = iV£;app(T), 



Eapp{T) — E - 



SE 



(15) 



which (via eqn. I14[) has a temperature independent stan- 
dard deviation equal to SE^/N indicating it to be a sta- 
tistically meaningful quantity for sufficiently large N. 



At the plastic transition temperature, the apparent 
barrier energy becomes 



-E-appCTc) — E 



1-2{SE/Efa. 



(16) 



Eqn. [in] is plotted in fig. [5]d as a function of a for the 
special case of £; = SE. The decrease of E ' ' (Tc) with 
respect to a reflects that Tc decreases with a. 

As a — >■ 0, the apparent barrier energy, i?app(7c), lim- 
its to NE which is the most probable barrier energy of 
P{E). On the other hand, when a equals one half, the 
apparent barrier energy becomes zero. These results can 
be understood by first realising that because of the ther- 
mal exponential factor in eqn. [S] low barrier energies will 
dominate the average from which i?app(?c) is obtained. 
For example, given M sampled barriers, it will be the 
lowest barrier energy of this list that makes the largest 
contribution to the apparent barrier energy. Thus the 
statistics of extreme values naturally emerge: sampling 
the distribution once (with a ~ 0) results on average 
in the most probable value of the distribution (which is 
NE), whereas sampling the distribution an increasing 
number of times will lead to lower barrier energies dom- 
inating the apparent barrier energy. The fact that a is 
limited to being less than or equal to one half is an ar- 
tifact of the Gaussian distribution not limiting to zero 
as the barrier energy limits to zero, giving £'app(2c) = 
at S = 1/2. A distribution for which P{E 0) ^ 
would allow for an infinite range of a and a finite posi- 
tive i?app(Tc) resulting in an everlasting approach to the 
zero barrier energy with increasing a. In sec. Ill C[ a mod- 
ification to the Gaussian distribution will be investigated 
that gives this correct limiting behaviour. 

Eqn. [T3l suggests that eqn. [9] may be written as 



M ( exp - 



-Mapp(T) exp - 



NE.ppjT) 



ksT J / '^'^^ ' knT 

where Afapp(T) is given by 

Mapp(r) = exp(iVaapp(r)). 



with 



AT) 



SE 
kWf 



(17) 
(18) 

(19) 



In eqn. 1171 the first factor is interpreted as the appar- 
ent number of structural transformations and the second 
factor is associated with the apparent barrier energy of 
eqn. 1151 At high enough temperatures, the apparent en- 
ergy limits to NE, the most probable value of the barrier 
energy distribution. Thus in the high temperature limit, 
the entire distribution contributes to the barrier kinet- 
ics, a fact that is refiected by the apparent number of 
available structural transformations limiting to M — the 
total number of structural transformations within each 
heterogeneous volume element. As the temperature low- 
ers, both the apparent barrier energy and the apparent 
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number of available structural transformations reduce be- 
cause, with a lower temperature, the higher barrier ener- 
gies become increasingly unlikely to occur. Eqn.fTOl shows 
that if T > Tc then eqn. [T7] will be exponentially large 
because the thermal factor is over-compensated by the 
apparent number of structural transformations. When 
T < Tc this turns to an under-compensation, resulting in 
eqn. 1171 becoming neglible. 

Is the functional form of eqn. [9] and [IT] valid at all 
temperatures below Tc? 

Inspection of eqns. [18] and [TU reveals that at a low 
enough temperature the argument within the exponen- 
tial will equal zero resulting in the apparent number of 
available structural transformations equalling unity. This 
will occur at the temperature 



Tf = 



SE 1 



(20) 



When T < T{, the apparent number of available struc- 
tural transformations will be less than one. This reflects 
that at this temperature some heterogeneous volumes 
will contain no available structural transformations with 
a barrier energy equal to E^^ppiT), and the dominant bar- 
rier energy will (on average) be the lowest available one, 
£'app(T'f), which equals 



(21) 



E,pp{Tf) ^ E - 5EV2a 



Upon further decrease in temperature, iSappCTf) will con- 
tinue to dominate where now the correct expression for 
T < Tf is no longer eqn. [9] or [TT] but rather 



M (exp [ — 



= cxp 



NE, ppjTf 
knT 



(22) 



Both Tf and i?app(Tf) are plotted as a function of a 
in fig. [2] showing that Tf and i?app(Tf ) generally occur 
below their plastic transition counterparts. 

Analogous low temperature behaviour is also seen in 
spin-glass systems, however with respect to thermody- 
namic variables where at a sufficiently low temperature, 
Tf , the thermodynamic entropy of the spin-glass becomes 
zero resulting in a temperature independent free energy 
upon further decrease in temperature. This is referred to 
as freezing, where at T < Tf one configuration of the sys- 
tem dominates the statistics giving an entropy of zero. 
This was first derived for the Random Energy Model 
by Derrida&8 using a micro-canonical ensemble approach. 
Because, eons. lUl and [T71 have the structure of an average 
partition function, the mathematical formalism associ- 
ated with the micro-canonical approach of Derrida can be 
also used as an alternative derivation of the right-hand- 
side of ean.[T71 and is presented in appendix [XI Carrying 
out this analogy allows the apparent barrier energy to 
be viewed as an internal barrier energy and the natural 
logarithm of the apparent number of available structural 
transformation as a barrier entropy, resulting in eqn. 1171 





temperature 



temperature 



FIG. 3. Plot of a) free barrier energy and b) plastic rate as . 
function of temperature. 



being written as 



M { exp 



E 
k^ 



cxp 



FiT) \ 
kBT J 




(23) 



(24) 



is the free barrier energy. For T < Tf , F{T) = T(Tf ). 

Within this analogy, the plastic transition temperature 
Tc represents the temperature at which the free barrier 
energy is zero — a condition that reflects equal contri- 
butions from the internal barrier energy and the barrier 
entropy. In other words, Tc corresponds to the transition 
from an internal barrier energy dominated to a barrier 
entropy dominated plastic rate^^ . It must be emphasised 
that the internal barrier energy, barrier entropy and free 
barrier energy are not true thermodynamic variables, but 
rather represent an alternative derivation and interpre- 
tation of eqn. 1171 Such an analogy also gives a proper 
context to the meaning of Tf , eqn. I20[ which will be now 
referred to as the kinetic freezing temperature — a tem- 
perature below which the barrier statistics are dominated 
by a single barrier energy scale. 

Figs. [3^ and b plot the free barrier energy, eqn. [^ 
and eqn. [22] (or eguivalentlv [T7)) . taking eqn. [52] into ac- 
count, as a function of temperature where both the Tf 
and Tc are indicated. In flg. ^jp, the proper treatment of 
the low temperature behaviour becomes apparent, where 
with decreasing temperature eqn. 1231 becomes a constant 
when T < Tf. Without such a correction, decreasing the 
temperature below Tf would result in the unphysical in- 
crease of eqn. [53] and therefore the plastic rate. Fig. [51d 
also demonstrates that the sharpness of the transitions is 
controlled by N. 



C. The application of a load 

The application of an external stress will affect the 
distribution of barrier energies via a change in its first 
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and second moments, {E) and {E'^). What should the 
stress dependence of these moments be? For the present 
work, only a pure shear stress will be considered. 

The occurrence of a structural transformation will re- 
sult in a localized rearrangement of atoms, which in the 
far-field limit can be accommodated by a homogeneous 
distortion of the materiali^i^. Thus the irreversible 
structural transformation will change the global stress 
state of the system. Since the zero-load P{E) has no 
knowledge of an external load, any structural transforma- 
tion corresponding to barrier energy E could with equal 
chance aid or hinder the applied stress state. In a ther- 
mal activation picture of plasticity, this corresponds to 
the equal chance of a particular barrier energy either in- 
creasing or decreasing upon the application of the pure 
stress. Hence, a broadening of the distribution of barrier 
energies will occur. The present work therefore consid- 
ers a stress independent first cumulant and a shear stress 
dependent second cumulant. The special case of a stress 
dependent first and second cumulant such that E = SE is 
maintained was considered in ref4^ and results in a yield 
criterion that is operationally similar to that derived by 
Johnson and Samwer'*^. 

By symmetry, the broadening should be an even func- 
tion of the pure shear stress, giving the leading order 
stress dependence of the second moment as^ 



(25) 



SE{a) = SE{0) I 1 + ( — 



with SE(0) equal to what was previously referred to as 
SE. The parameter ctq will be discussed in more detail 
in sees. IIIII and IIVI Use of eqn. [25] results in the yield 
transition temperature, eqn. 111! becoming 



Tcia) 



E 



{SE{0)/E)'{l + {a/aory 



1-^1-2 {SE{0)/eY (1 + {a/ao)^fa 



(26) 

and the freezing transition temperature, eqn. I20| becom- 
ing 

6E{0) (1 + i<j/<Jo?) 



Tt{a) 



2a 



(27) 



Thus both Tc and Tf become functions of the applied 
stress. 

Inspection of eqns. and [57] reveal there exists a 
stress, (Tf, at which the plastic transition temperature 
and the freezing temperature are equal. Indeed, T'c(o't) = 
Tf(crf) gives 



1 



/2S 



(28) 



and thus A:B'7c(crf) — kBT[{a{) = E/2a, with the appar- 
ent barrier energy at freezing equal to zero (eqn. [2T|). 

Fig. [3^ plots the applied shear stress, a in units of cto, 
versus the plastic transition temperature Tc (a) (eqn. 
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FIG. 4. a) Plot of plastic transition stress as a function tem- 
perature for values of a ranging between 0.1 and 0.4. The 
dashed lines represent the corresponding stress dependent 
freezing temperature, b) Shows the same curves scaled with 
respect to the kinetic freezing stress, erf and the plastic transi- 
tion temperature at zero- load, Tf . c) Shows the apparent and 
freezing (dashed lines) barrier energy and d) the apparent a 
which gives information on the number of thermally accessible 
structural transformations. Both c) and d) are with respect 
to temperature scaled with the plastic transition temperature 
at zero- load, Tf. 



and the freezing temperature Tf(a) (eqn. [27]) for a range 
of a. In fig. [4] the special case oi E = SE{0) is consid- 
ered. Fig. [3^ shows that by increasing the applied shear 
stress, Tc{a) reduces from its zero load value, whereas 
T[{a) increases, eventually resulting in Tc{(J[) = T[{a[) 
at the kinetic freezing shear stress a^. Below this tem- 
perature eqn. [5S] is no longer applicable and the plas- 
tic transition rate becomes a constant given by eqn. 1221 
Upon increasing a both the maximum yield stress and 
zero-load plastic transition temperature reduce. When 
plotted as scaled quantities, fig. the plastic transi- 
tion temperature curves close to T/Tc depend relatively 
weakly on the choice of a. Fig. [4]; shows the apparent 
barrier energy and kinetic freezing barrier energy at the 
plastic transition stress as a function of the plastic tran- 
sition temperature. It is seen that the apparent barrier 
energy is at a maximum at zero load and then reduces 
in an approximately linear fashion with reducing tem- 
perature, finally reaching zero at freezing. The freezing 
barrier energy also reduces (non-linearly) to zero at freez- 
ing. Fig. [4ji plots aapp{T) as a function of the plastic 
transition temperature, showing a non-linear reduction 
with respect to decreasing temperature, to zero at freez- 
ing, demonstrating that the apparent number of available 
structural transformations rapidly reduces to zero as the 
temperature decreases from Tc. 
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The picture which therefore emerges, is that at zero 
load, there exists a temperature where a sufheient num- 
ber of transitions exist to offset the very small thermal 
factor corresponding to the apparent barrier energy at 
that temperature. This gives a non-negligible plastic 
rate. At a lower temperature, upon the application of 
a high enough applied shear stress, such a non-negligible 
plastic rate can again be achieved. However, since the 
temperature is lower, both i?app(?c) and aapp(T'c) have 
reduced due to the higher barrier energies of the distri- 
bution becoming thermally inaccessible. Due to this re- 
duction in the available number of structural transforma- 
tions, there exists a temperature Tc(crf) — T[{a[) below 
which the system kinetically freezes before the plastic 
transition is obtained, irrespective of the magnitude of 
the applied shear stress. By increasing a (and therefore 
increasing the total number structural transformations), 
both the temperature and stress at which non-negligible 
plasticity and/or kinetic freezing is achieved, is reduced. 

Thus, with the use of a Gaussian distribution, some 
features of the high temperature regime of deformation 
of BMGs are qualitatively obtained. Fig. HJd shows that 
the plastic transition stress rises rapidly from zero when 
the temperature is decreased from the zero-load plastic 
transition temperature, and then saturates at a transi- 
tion temperature coinciding with the applicability limit 
of the Gaussian model. In what follows it will be shown 
that, when a more realistic distribution of barrier energies 
is used, the developed theory may be extended to lower 
temperatures with a corresponding change in deforma- 
tion behaviour that may be associated with the experi- 
mentally observed low temperature deformation regime 
of BMGs. 



D. A more realistic distribution of barrier energies 

That the freezing regime is entered in fig. 2] is a direct 
result of the Gaussian distribution having a finite (but 
small) probability at the zero barrier energy. Because 
of this, both the apparent barrier energy and freezing 
barrier energy become stuck at a zero barrier energy. To 
modify the distribution in a way that has the correct limit 
P(i? — >■ 0) — >■ and still facilitates the developed math- 
ematical formalism, the argument of the exponential of 
the Gaussian distribution must be modified. 

Starting from 



P{E) - exp 



N /E/N - 

T I W 



E 



such a modification would entail 
P{E) - exp 



TV f g{E)/N-E 
2 \ 6E 



(29) 



(30) 



where the function g{E) is chosen such that the resulting 
apparent and freezing barrier energies remain extensive 
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FIG. 5. Similar data to that of fig. [4] using the modified 
Gaussian distribution of barrier energies where both a high 
and low temperature regime of plasticity become apparent. 



quantities. Additionally, g{E — 0) — >■ cxd so that P{E —>■ 
0) ^ 0, and at larger values of E, g{E) c:^ E to retain 
the Gaussian form when E is comparable to E. A simple 
choice that satisfies these requirements is 



giE) = E- 



E 



(31) 



where a is a parameter with units of energy. 

Insight into the consequences of such a modification 
and into the meaning of a may be obtained by writing 
the resulting stress dependent apparent barrier energy at 
freezing as (see appendix \^ eqn. lAlSp 



Etia) 



Ef{a) 



(Efia))' 



4a2 



(32) 



Here Ef (a) is the stress dependent barrier energy at 
freezing when using a pure Gaussian distribution, with 

E{ (dp) ~ 0. The above equation reveals that with the 
modified Gaussian, the barrier energy at freezing will 
never become zero and at the stress, crp, at which freezing 
occurs for the pure Gaussian model, Ef{a^) = a. Thus 
the parameter a sets the apparent barrier energy at the 
freezing stress originally given by the Gaussian model. 

To evaluate the integral of eqn. [9] for this modified 
Gaussian, the procedure outlined in appendix |A] is used 
to construct the so-called free barrier energy, F{T) from 
which, F(Tc) = can be numerically solved to obtain the 
critical temperature at which yield occurs for zero load 
conditions. 

Fig. [5] displays similar data as that shown in fig. 01 In 
this figure a value of a = 0.05i? was used. Inspection of 
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fig. [5] reveals that at a given temperature, kinetic freez- 
ing will always occur at a stress lower than that of the 
plastic transition. Thus a plastic transition temperature 
will exist for all temperatures below the zero load plastic 
transition temperature. Moreover, use of the modified 
Gaussian distribution gives three distinct temperature 
regimes of behaviour: 1) a high temperature regime, in 
which the stress at which the plastic transition temper- 
ature occurs rises rapidly and is well described by the 
pure Gaussian (see fig. |3]of sec. IIIGp . 2) a low temper- 
ature regime in which the plastic transition stress rises 
much more weakly (approximately linearly) with decreas- 
ing temperature and 3) an ultra-low temperature regime 
in which the plastic transition stress diverges. Regime 3) 
occurs at low enough temperatures and high enough plas- 
tic transition stresses, such that the thermal activation 
hypothesis becomes questionable and where an athermal 
regime of plasticity would emerge, a phenomenon that is 
beyond the scope of the present work. For the low tem- 
perature regime, 2), figs.[SJ: and d, demonstrate that the 
modified Gaussian averts a zero apparent barrier energy 
and aiipp{T) at finite temperature, giving a distinctly dif- 
ferent temperature dependence to that of the high tem- 
perature regime. Decreasing a is found to reduce the 
gradient of the approximately linear temperature depen- 
dence of the plastic transition stress, where in the limit 
of small a a plateau is reached for temperatures less than 
the pure Gaussian derived Tc(crf) — Tf (af) value. 

In summary, with the modification, P{E — 0) — > 0, 
the present thermal activation model is now able to cap- 
ture two distinct temperature regimes of deformation — 
with increasing temperature the low temperature regime 
exhibits a weak decrease of the plastic transition stress, 
with a smooth and analytical transition to the high tem- 
perature regime of deformation where the plastic transi- 
tion stress drops rapidly to zero at a finite temperature. 
In the next section, it will be demonstrated that the de- 
veloped theory can quantitatively reproduce the corre- 
sponding experimental trends of some well known BMG 
materials. 



where now eqn. [7]is used. Doing so, gives 



III. 



APPLICATION TO THE EXPERIMENTAL 
REGIME 



In sees. Ill Bi to III Dl the plastic transition stress was de- 
rived from the condition M {exp{— E / (k-BT))) — 1. For 
application to a real glass material, the full equation for 
the characteristic plastic rate (eqnl?]) must be used. The 
plastic transition stress must also take into account the 
timescale, Texp, of the deformation experiment. In con- 
sidering the experimental timescale, a more appropriate 
definition of non-negligible flow to occur would be when 
the timescale associated with plastic activity is equal to 

'^Gxp : 



1 

N 



In 



^pOO 



pOO 



knTr 



1 / SE 



E 



kuTr 2 V kuT, 



I f 5E '-^ 



E 



knTr, 2 V k-aT, 



-0, 



with 



and 



llnf^") 
N V ^poo J 



e' = E-' 



N 



^ 
(34) 

(35) 
(36) 
(37) 



M"'(Tc)^l, 



(33) 



Therefore the inclusion of an experimental timescale and 
the prefactor, Tpo(T), into the plastic transition condi- 
tion renormalizes the parameter a and the mean barrier 
energy E. For kinetic freezing, the corresponding bare 
parameters should be used. 

To apply the present model to a real material, it is 
recognised from figs. |4] and [5] that the high temperature 
regime is well described by a pure Gaussian barrier en- 
ergy distribution. Since simple analytical formulae for 
all relevant quantities are possible for a pure Gaussian, 
the high temperature deformation properties will be used 
to give estimates of all but one model parameter. The 
model parameters to consider are 

1. a, the configurational barrier entropy which gives 
the total number of barriers per heterogeneous vol- 
ume. 

2. N, the number of atoms per heterogeneous volume. 

3. E and SE, the mean and standard deviation per 
atom of the barrier energy distribution. 

4. the logarithm of Toxp/Tpoo and -Epoo- 

5. (To, the shear stress sensitivity of the distribution 
broadening. 

Apart from Toxp , all of the above parameters are related 
to the material properties of the particular structural 
glass of interest. 

Inspection of figs. |4] and [5] re veal three quantities which 
define the range of the high temperature regime: the 
plastic transition temperature at zero load, and the tem- 
perature and stress at which the low temperature regime 
is entered (for the pure Gaussian when Tc{a[) — T[{a[)). 
These three quantities may be determined directly from 
experiment. Indeed for a sufficiently long enough experi- 
mental timescale, Texp, the plastic transition temperature 
at zero load will be close to the material's glass transition 
temperature. This fact will be exploited in what follows. 
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The temperature at which the material enters the low 
temperature regime of plasticity is also experimentally 
well defined, being characterised by the transition away 
from thermally activated viscoplasticity which usually 
occurs at a temperature CS-O-OTg = 7Tg = Tc{ai) — 
T{{a{), and a plastic transition stress of erf. Finally, the 
sharpness of the plastic transition (set, in part, by TV) is 
also an accessible experimental quantity via the kinetic 
fragility^, which is defined as 



^Kinetic 



d log 77 



d{T,/T) 



(38) 



T=r„ 



Here the viscosity, ij may be written via the Maxwell re- 
lation as 77 = GooTp{T), where G is the instantaneous 
shear modulus and Tp (T) is now viewed as the relaxation 
time scale of the under-cooled liquid. The kinetic fragility 
gives a dimensionless measure of how rapidly the viscos- 
ity drops upon a small increase of temperature from the 
glass transition temperature. At this temperature one 
can equally well speak of a viscosity whose associated 
time scale is the inverse of the plastic rate. Doing so 
results in a simple formula for the kinetic fragility 



TTT'kinctic 



^^app(Tg) 

khT„ 



(39) 



where again Tc ~ Tg is assumed. Hence the kinetic 
fragility sets the apparent barrier energy of the het- 
erogeneous volume at the glass transition temperature. 
Eqn. [551 is ^ well known result, when viewing fragility as 
a kinetic phenomenon, however that it emerges from a 
time scale derived from the product of two exponentials 
whose arguments both depend on temperature (eqn. [T7)) 
is not obvious and constitutes a clear justification of in- 
terpreting i?app(T') as the apparent barrier (or activation) 
energy. _ 

\ of the above allow E, SE{0) and N to be deter- 
mined as a function of a, E'poo and In (rcxp/rpoo). Indeed 
by exploiting Tc ~ Tg, and eans.E51and[5^ the following 
relations may be obtained 



NE = 



NdE{0) = 



Na 



1-7 

27-1 
1-7 

1 

2(1-7) L 



mloglO-ln(^^^"P 



, TpOO 



^loglO-^— Inf^^ 
1-7 V ^POO / 



(40) 



mlogl0-(2-7)lnf^) 



ctq is given by 



Of 



y/2Na"/ i?^^ 1 

y/NSE{0) 



(41) 
(42) 

(43) 



where typically, ai, is equal to 0.02G with G being a 
representative (pure shear) elastic modulus. 
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FIG. 6. a) Plastic transition and kinetic freezing stresses as a 
function of temperature, b) the apparent and kinetic freezing 
barrier energy per homogeneous volume as a function of tem- 
perature, c) Apparent number of available structural trans- 
formations as a function of temperature and d) the plastic 
transition rate as a function of applied shear stress for the 
temperatures 300K and 610K. 



The above approach is motivated by the known strong 
correlation between the characteristic barrier energy 
scale of the a relaxation mode (derived by applying an 
Arrhenius viscosity law to eqn.[38]and using experimental 
values for Tg and mkinetic) and the activation energy of a 
simple thermal activation model of plasticity that is able 
to describe well the high temperature (T > O.STg) de- 
formation properties of glassesi^ii^i^. The above is also 
consistent with the initial assertion that the mega-basin 
barrier energy landscape is associated with the a relax- 
ation mode. Furthermore, by associating EpQQ with the 
/3 relaxation mode, this material parameter may also be 
directly determined from experimental measurements of 
the /? barrier energy via differential-scanning-calorimetry 
or dynamical-mechanical-spectroscopy methods^^— . 

The remaining parameter is the ratio Toxp/Tpoo. T'oxp 
is related to a (sufficiently low) reference strain-rate, e. 
To first order, their relationship will be texp oc [i] , 
where the proportionality constant will depend on the 
details of the characteristic strain associated with a meta- 
basin escape and on temperature. These aspects, vis d 
vis a stress-strain relation, will not be considered in detail 
in the present work, however appendix |B] gives a simple 
estimate of the proportionality constant resulting in the 
simplified relation texp — 0.01 [e] ^. Since, the escape 
from a mega-basin is assumed to be mediated by /3 mode 
relaxation, the Tpoo will be some multiple of the system's 
Debye frequency. Presently it will be taken as Tpoo = 
1 X 10~^"'^/sec. Fortunately, it is the logarithm of the 
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ratio of Tcxp/Tpoo that enters into the model, making the 
model quite insensitive to the precise orders of magnitude 
of these timescales. 



Parameter 


Taken experimental value 




623 K 


^^kinctic 


50 


i;"(=ln(10)mkincticTg) 


6.1 eV 




1.43 eV 


Shear modulus 


34.1 GPa 



TABLE I. Parameters taken from refs.— of Vitreloy-1 used 
to determine the model parameters. 

Fig. shows the yield stress versus temperature for 
the experimental quantities of Vitreloy-1 (see tab. |I]) 
at a strain rate of £ = 10~^/sec, which is assumed to 
give a zero-load transition temperature close to the glass 
transition temperature. In these curves, 7 = 0.9 and 
Na ~ 0.125 eV. Inspection of this figure shows that 
the plastic transition stress rises rapidly with decreas- 
ing temperature in the regime T > O.OTg, whereas for 
T < O.QTg a change in behaviour occurs resulting in 
a weaker, approximately linear increase in yield stress 
with decreasing temperature. These curves show that at 
the plastic transition stress, the phenomenon of freez- 
ing is clearly avoided. Fig. ^jp displays the apparent 
barrier energy as a function of temperature. At the 
zero- load limit, the apparent barrier energy has a value 
that is equal to the experimental value of E" — E^^ = 
mki„eticfcBTgln(10) - Epoo - 6.1 - 1.43 eV. With de- 
creasing temperature/increasing yield stress, the appar- 
ent barrier energy rapidly reduces until the low tem- 
perature regime is reached. In this regime, the plastic 
transition stress increases approximately linearly with re- 
spect to decreasing temperature. Fig. [BJa also displays 
the freezing barrier energy. Both barrier energies remain 
positive for all temperatures. Fig. [SJ; shows the natural 
logarithm of the apparent number of available structural 
transformations, aapp{T)N and, as with the apparent 
barrier energy, its value rapidly reduces and saturates as 
the low temperature regime is entered. In this regime, 
its value changes little. 

In the low temperature regime of deformation, experi- 
ment also reveals a robust elastic regime of stresses and a 
sharp transition to plasticity — behaviour which is also 
seen in the present model. Fig. [6ji plots the plastic rate 
as a function of applied shear stress for the two temper- 
atures 300K and 61 OK, where the latter is clearly in the 
high temperature regime of deformation where no kinetic 
freezing transition stress regime is encountered, and there 
is a more gradual transition to plasticity. However at a 
temperature of BOOK, much of the stress regime below the 
plastic transition stress occurs within the kinetic freez- 
ing regime and the transition to plasticity is considerably 
sharper. 

All curves are found to be independent of a. What is 
the origin of this independence? Inspection of eqns. HD] 
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FIG. 7. a) Plot of N, the characteristic number of atoms 
per heterogeneous volume and b) the mean and standard de- 
viation per atom of the barrier distribution. Both figures 
are derived from the parameters of Vitreloy-1, tab.]!] c) Bar- 
rier energy distribution for the unloaded state and at a ap- 
plied stress of 0.02G, for the parameter set corresponding to 
Viterol-1. 

to [151 reveals the quantities on the left hand side are 
fully defined via the experimental parameters of the right 
hand side. Thus the mean NE (eqn. |40|) and standard 
deviation ^/NSE (eqn. [41]), and thus the Gaussian part 
of the distribution, are fixed by measurable deformation 
properties and independent of a. On the other hand, 
choosing an actual value of a will fix N via eqn. [12] and 
therefore determine E and SE. 

Thus the experimental temperature dependence of the 
plastic transition stress uniquely determines the underly- 
ing distribution of a-relaxation barrier energies, whereas 
a sets the value of N, and therefore the size of the het- 
erogeneous volume. This must be viewed as a setting 
N, such that M = exp(aV) is of a sufficient size for the 
observed phenomenon to occur. It is in this way that the 
deformation properties are insensitive to the underlying 
microscopic detail defined via a. Fig.[7K plots eqn.l42land 
shows that for values of a < 1, N can be of the order of 
several thousand atoms, giving corresponding values for 
E and SE in the tens of meV per atom (fig. [TId). This 
gives a mean barrier energy of NE — 36.5 eV and stan- 
dard deviation ^/NSE = 1.3 eV. Fig. [7]; plots the result- 
ing a-relaxation barrier energy distribution, also show- 
ing the regime of barriers accessed for T < Tg, as seen 
in fig. [BJd. Inspection of this figure demonstrates that 
the corresponding deformation properties probe only the 
extreme low barrier energy tail of the distribution, with 
the most probable part of the distribution (giving the 
dominant contribution to M) playing little role. Thus 
in the low temperature regime, the thermally accessi- 
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FIG. 8. Replot the plastic transition stress as a function of 
temperature in fig. [S^, where temperature is scaled with re- 
spect to the glass transition temperature. Experimental data 
for Viteroly-1 and other bulk metalic glasses are also shown 
(blue data taken from Johnson and Samwer*^ and red data 
from Lu et af-). 

ble a-relaxation PEL has significantly flattened, with the 
coarse grained barrier energy scale being comparable to 
that of the characteristic /3-relaxation/STZ barrier en- 
ergy. 

Fig. [5] displays the typical experimental yield stress 
versus temperature data used to determine the optimal 
value of a. In this figure, the data from both experiment 
and the present theory is plotted as a yield stress (scaled 
by a representative shear modulus) versus temperature 
(scaled by the glass transition temperature). Shown is 
experimental data for Vitreloy-1 and other bulk metallic 
glasses, taken from Johnson and Samwer^l (presumably) 
at a number of different strain rates. Also displayed is 
Vitreloy-1 data at a strain rate of 1 x 10~^/sec spanning 
both the low and high temperature regimes of deforma- 
tion^. The figure shows good agreement between the 
current model and that of the Vitreloy-1 experimental 
data, and also emphasises the experimentally observed 
universal temperature dependence of the low tempera- 
ture deformation regime. Due to the overly simple model 
for plastic strain (appendix IB|) used to obtain a relation 
between the characteristic experimental time scale and 
the strain rate, the strain rate dependence of the devel- 
oped model will not be considered in the present work. 

IV. DISCUSSION 

The present model has some similarities with the well 
known thermodynamic model of Adam and Gibbs^ for 
undercooled liquids. Both assume an exponential number 
of available structural transformations, although in the 
case of Adam and Gibbs, it is assumed that this num- 
ber is equal to the number of accessible final states of 



the system, allowing Sconfig{T) in exp (sconfig(7')A''*) to 
be associated with the configurational entropy per parti- 
cle. Here N* is the size of the cooperatively rearranging 
regions (CCR) whose temperature dependent value is set 
via the assumption that Sconfig(2^)Af* is equal to a given 
fixed value. By assuming that the activation energy is 
proportional to the size of the CCR, and therefore N*, 
the activation energy becomes inversely proportional to 
the configurational entropy per particle leading upon fur- 
ther assumptions to the Vogel-Tammann-Fulcher equa- 
tion^"— for the characteristic relaxation time of an un- 
dercooled liquid. 

For the presently developed model, the analogous a is 
a fixed material constant and is referred to as the con- 
figurational barrier entropy per atom, N is the charac- 
teristic number of atoms within a heterogeneous volume 
element, and the product aN is set by the overall exper- 
imental temperature dependence of the plastic transition 
(yield) stress close to the glass transition (eqn. 142]) as are 
the first and second cumulants of the barrier energy dis- 
tribution (egns. HUl and HT|) . Here, the apparent barrier 
energy is derived from the thermally accessible part of 
the barrier energy distribution, and although this barrier 
energy is also dependent on N, its temperature depen- 
dence is quite different from Adam and Gibbs, decreasing 
rather than increasing as the temperature drops. Indeed, 
in the present work, the rapid increase in the relaxation 
timescale (the inverse of the plastic transition rate) as 
the temperature decreases arises from the strong reduc- 
tion in the apparent number of states and therefore clapp- 
From this perspective a more useful comparison might be 
between Sconfig(2^) and aa_pp{T), since both decrease with 
temperature and both arc ultimately responsible for the 
rapid increase in the relaxation timescale as the temper- 
ature reduces. 

The present model may be partly viewed as a deriva- 
tion of a temperature dependent critical barrier energy 
given by -Epoo + NEa,pp(a,T{) which has contributions 
arising from both the a and /3 (STZ) relaxation modes. 
Such a critical barrier energy has been assumed by both 
Argon'^ and Johnson and Samwer— , although with no 
explicit temperature dependence, allowing for the use of 
a single barrier energy. Indeed for the present model, 
the derived temperature dependence is weak when com- 
pared to that of the high temperature deformation regime 
(figs. [St and [B]:). Together eqns. [T3] and [T3] (the first 
and second moments of the distribution P{E) exp{—/3E)) 
show that the a-relaxation mode contribution to such a 
critical barrier energy has standard deviation that scales 
as \/N, where N is the characteristic number of atoms as- 
sociated with the heterogeneous volume-scale of a BMG. 
For small enough a, N can be of the order of a few 
thousand, resulting in N Eiipp{a,T{) being a statistically 
meaningful quantity and therefore justifying the notion 
of a single critical barrier energy. Although eqns. [13] and 
1141 are formally only applicable to the high temperature 
regime where a simple Gaussian suffices, numerically it 
was found that this is also the case for the low tem- 
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perature regime of deformation when using the modified 
Gaussian distribution. 

The current work, shows that in addition to a criti- 
cal (apparent) barrier energy, there also exists a critical 
(apparent) number of available structural transformation 
and together these probe the extreme low energy tail of 
the barrier energy distribution (fig. [71:). When T > Tf, 
this extreme low energy tail is characterised directly by 
the thermal distribution, P{E) exp{—(3E), however for 
temperatures at and below the kinetic freezing tempera- 
ture, Tf, the statistics changes. In this regime, one bar- 
rier energy scale dominates and is equal to NEapp{Tf). 
This kinetic freezing barrier energy is an average value, 
and when sampling heterogeneous volumes a distribu- 
tion of such kinetic freezing barrier energies will be ob- 
tained. For positive valued distributions such as P{E) 
the Fisher- Tippett-Gnedenko (FTG) theoreniZ2.i23 states 
that such an extreme value distribution will be of the 
Weibull form''^. Such a theorem for extremal values is 
analogous to the well-known central limit theorem for 
averages, and does not depend on the precise form of 
P{E). The connection between the phenomenon of freez- 
ing and extreme value statistics has been established 
by Bouchard and Mezard^^, demonstrating that like the 
statistics of the extreme, the statistics of freezing be- 
longs to a particular universality class suggesting a ro- 
bustness against material specifics. Appendix [C] details 
a similar connection between the present kinetic freezing 
phenomenon and extreme value statistics, showing that 
the freezing barrier energy derived in sec. |TT] is actually 
an underestimate of the extreme value statistics derived 
freezing barrier energy — their difference arising from 
the finite size of the heterogeneous volume. 

The two distinct temperature regimes of deforma- 
tion exhibit clear differences in temperature dependence. 
What is the nature of their difference in terms of the 
underlying plastic activity? In the high temperature 
regime, each heterogeneous volume will contain a suffi- 
cient number of thermally accessible structural transfor- 
mations such that plastic activity always occurs — the 
degree of plasticity over the experimental time-scale Toxp 
of any sampled heterogeneous volume is well reflected by 
the derived average. In this regime, all heterogeneous 
volumes are deforming resulting in the homogeneous on- 
set of plasticity. 

To gain insight into how the low temperature regime 
differs from the above scenario, the kinetic freezing limit 
is first considered. At a low temperature, T, and at an 
applied stress, cr, close to the freezing stress, the average 
plastic rate is given by 

Tp ~ Tpoo X exp ( ^- 1 , (44) 

which represents the effect of only one available struc- 
tural transformation within the heterogeneous volume. 
In this regime, sampling a particular heterogeneous vol- 
ume over the time period Texp can result (with a par- 
ticular probability) in that volume not deforming. Thus 



the (negligible) plasticity at the kinetic freezing limit be- 
comes strongly heterogeneous when sampling heteroge- 
neous volumes. Raising the applied stress exits the ki- 
netic freezing regime and the plastic rate is now given 
as 

Tp ~ Tpoo X Afapp(CT, T) X exp ( ^ ) ' 

(45) 

Inspection of figs. [6}d and c reveal that at room tem- 
perature for Vitreloy-1, Mapp((T, T) is of the order of 100, 
whereas i?app(o', T) has not increased greatly from that 
of its freezing value. Although plasticity is more likely in 
any sampled heterogeneous volume at the actual plastic 
transition stress, in transiting to non-negligible plastic- 
ity, the statistics is expected to be strongly influenced by 
those of the kinetic freezing regime. 

Thus the transition from the high temperature regime 
to the low temperature regime is characterised by a 
change of statistics from that where the average plas- 
tic time-scale reflects well the degree of plasticity of 
any sampled heterogeneous volume to that where the 
statistics of small numbers and extremal values admit 
a non-negligible probability of any particular heteroge- 
neous volume not deforming. The implication is that 
in passing from the high to low temperature regime, the 
plasticity becomes inherently inhomogeneous in the tran- 
sition from elasticity to plasticity — the so-called micro- 
plastic regime of deformation. How such a low tempera- 
ture heterogeneity manifests itself as an emerging mate- 
rial instability in the macro-plastic/flow via (say) shear 
banding would give fundamental insight into the transi- 
tion from homogeneous to heterogeneous plasticity seen 
experimentally as the temperature decreases from the 
glass transition. 

The generality of the current result is now discussed. 
Whilst the use of a modified Gaussian is certainly an 
approximation to the true distribution of a-relaxation 
barrier energies, the existence of both a kinetic freezing 
and plastic transition temperature is insensitive to the 
precise form of P{E), where its essential properties are 
the extensive first and second cumulants — a require- 
ment that is intimately related to the known exponential 
scaling of the number of stationary configurations of a 
structural glass. Moreover, at and close to the freezing 
regime, where extreme value statistics dominate, the re- 
sulting Weibull distribution via the FTG theorem is a 
universal result independent of the precise form of P{E). 
One implicit approximation thus far not discussed, is that 
the average over heterogeneous volumes is taken with the 
assumption that each such volume is statistically inde- 
pendent of the other, that is, they do not interact. A 
similar assumption is made for the CCR of the Adam 
and Gibbs formalism^. The assumption of a lack of 
an explicit interaction between heterogeneous volumes is 
clearly an approximation to make the problem analyti- 
cally simple. The universal phenomenon of freezing has 
however been found to be robust against certain forms of 
interactionsIS-iH, and it is expected that for the low tem- 
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perature regime where kinetic freezing plays a defining 
role in determining the nature of plasticity, the current 
results will not be fundamentally changed by the inclu- 
sion of interactions between the heterogeneous volumes. 
It remains a topic of future work to, in the first instance, 
characterise the interaction between heterogeneous vol- 
ume elements and then to implement numerical proce- 
dures that are able to investigate their role not only in 
the transition from elasticity to plasticity but also in the 
flow regime of macro-plasticity. 

The Arrhenius form used for Tpo, eqn.[6l to describe the 
average properties of the atomic scale /3-relaxation/STZ 
processes mediating the a-relaxation landscape explo- 
ration is the simplest choice that embodies the thermal 
activation hypothesis. More complex temperature de- 
pendencies could be envisaged emerging from the distri- 
bution of barrier energies known to also exist for STZ 
barrier energies derived from atomic scale PEL explo- 
rationaiSiiiiZ^ using the ARTn technique^^i"— . Due to 
the independence of such phenomena to the size of the 
heterogeneous volume, N, sec. IIIII shows however that 
any such temperature dependence will only renormalize 
the parameters of the current model to a leading order 
of 1/A^. Thus the present model should be insensitive 
to the precise form of Tpo whether it be a function of 
temperature and/or applied stress. 

The current work has only considered an applied stress 
that is a simple shear. In comparing to experiment, it 
was found that the a representative shear stress modu- 
lus sets the proportionality constant between the applied 
shear stress and the broadening of P{E), eqn.H51 Defor- 
mation experiments are however usually performed under 
uni-axial loading conditions, were a compressive/dilatory 
component is also present. A more complex applied stress 
that contains such an isotropic component, such as the 
uni-axial loading experiment is expected to affect P{E) 
in a more complex manner by depending not only on a 
pure shear modulus but also on the Poisson ratio. This 
aspect, and thus the role of local volume changes, will be 
investigated in future work. 



V. CONCLUDING REMARKS 

In the 2000 review article of Angell et afi^ , which artic- 
ulated the contemporary questions to be addressed in the 
field of structural glasses, one such question for the low 
temperature regime was "What types of processes remain 
active in the glass when the a relaxation has been com- 
pletely frozen?" . The present work demonstrates how 
such a freezing of the a relaxation potential energy land- 
scape could occur as the temperature reduces. Indeed as 
the temperature decreases an increasing portion of this 
landscape becomes thermally inaccessible and therefore 
frozen — a process that is only complete at a temperature 
of absolute zero. Thus the statistics of the a relaxation 
processes, so integral to the under-cooled liquid regime, 
are also found to play a central role in determining those 



thermally active processes that lead to low temperature 
plasticity. 

In summary, a model has been developed to describe 
the temperature dependence of the transition from elas- 
ticity to plasticity, the micro-plastic regime, in a bulk 
metallic glass. Central to this model, is the assump- 
tion of thermally activated plasticity, and in particu- 
lar, that plasticity occurs via the thermal activation of 
a-relaxation processes that are themselves mediated by 
multiple thermally activated /3-relaxation/STZ activity. 
Further, it is assumed that the known exponentially num- 
ber of stationary points in a glassy potential energy land- 
scape for the under-cooled liquid regime, results in a dis- 
tribution of a-relaxation barrier energies that has first 
and second cumulants which are extensive with respect to 
the underlying heterogeneous volume-scale of the glass. 

Two distinct temperature regimes emerge. At high 
temperature, the shear stress at which significant plas- 
ticity occurs rises rapidly from its zero value close to 
the glass transition temperature, corresponding to a drop 
in the apparent barrier energy and apparent number of 
available structural transitions. In the case of the latter, 
this reduction is comparable to the many orders of magni- 
tude seen in the increase in the viscosity as the glass tran- 
sition temperature is approached from the under-cooled 
liquid regime. At a low enough temperature, the ap- 
parent number of structural transformations saturates, 
due to the eventual onset of kinetic freezing, and the low 
temperature regime is entered where the plastic tran- 
sition stress now increases approximately linearly with 
respect to decreasing temperature. In this deformation 
mode, the transition from elasticity to plasticity is sharp 
with respect to the applied stress, and for applied stresses 
characteristic of the elastic regime plastic activity is ki- 
netically frozen indicating a single dominant barrier en- 
ergy resulting in an underlying robustness of the elastic 
regime. 
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Appendix A: An analogy to thermodynamics and 
the phenomenon of freezing 

Further insight into eqn. [S] can be obtained by recog- 
nising that 

where (• • ■ ) is an average over heterogeneous volumes, 
has the structure of an environmentally averaged parti- 
tion function average. Thus the statistics of the current 
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barrier energy problem can be mapped to an equilibrium 
statistical mechanics framework, which then allows for 
application of Derrida's micro-canonical approach to the 
phenomenon of freezing^^. 

The average number of barrier energies between E and 
E + dE, {n{E)), is given by 



{n{E)) = M{E)dE = MP{E)dE, 



(A2) 



where dE must be small enough to ensure a well defined 
barrier energy but also large enough so that {^^{E)) is a 
smooth function of E (see refi^ for a related discussion). 
For a sufficiently large heterogeneous volume, fluctua- 
tions in ean. lA2l become small and {il{E)) becomes a sta- 
tistically meaningful quantity, allowing for a correspond- 
ing barrier energy to be defined via S{E) = fee ln(f2(_E)). 
For a Gaussian distribution with mean NE and standard 
deviation \/N6E this is 



'^^)- (A3) 



+ kB In 



,V2TiN6E 



From equilibrium thermodynamics, an effective temper- 
ature is formally defined via 



dS{E) 



dE E=U(T) 

and from eqn. IA3[ this gives. 



1 



U{T)^N\E ^^^^ 



kBT 



(A4) 



(A5) 



as the relation between internal barrier energy and effec- 
tive temperature. It is noted that the internal barrier en- 
ergy, U (T) from eqn. IA5I is equal to the apparent barrier 
energy of eqn.[T3l With the barrier entropy (eqn. lASj ) and 
the barrier internal energy (via eqn. IA5|) . the free barrier 
energy, F{T) = U{T) - S{T)T, can be constructed: 



= -kjiTN a - 



E 



-f— y 

2 \kBTj 



(A6) 



where in the last step, egn. lASI has been used to make the 
free barrier energy dependent on effective temperature. 
Since the average partition function is proportional to 
the logarithmic of the free barrier energy, eqn. lAll can be 
written as 



M ( exp 



E 
k^ 



exp 



F{T) \ 
kBT ) 



(A7) 



which is identical to eqn. O That eqn. [S] can be achieved 
in this manner, indicates that the definition of an effective 



temperature via eqn. IA4l formallv corresponds to the true 
temperature of the system. 

In eqn. IA6[ the second term of eqn. IA3I has been 
dropped under the assumption that the heterogeneous 
volume is large enough such that it contributes neglibly 
to both the barrier entropy and the free barrier energy. 
This highlights the fact that the present derivation and 
that of sec. |TT] are formally equivalent only in the limit 
N — !■ oo. Continuing in this limit, at a sufficiently low 
enough barrier energy {^(Ef)) — 1 gives 



s(E,,.«J.-i(^ 



which has the relevant solution 




Ef^NlE- SEV2S 



= 0, (A8) 



(A9) 



or in terms of temperature {Ef — U{Tf )) via eon. IA4 



Ti 



SE 1 
fee V2^' 



(AlO) 



Eqns.EiandETO] are identical to those derived for the 
freezing barrier energy and temperature in sec. IIIBI 

From the perspective of thermodynamics, at the tem- 
perature Tf the entropy is equal to zero and since upon 
further decrease in temperature the thermodynamic en- 
tropy cannot further reduce, the system freezes into this 
single dominant state resulting in a temperature indepen- 
dent free energy^. For the present analogy to barrier en- 
ergy kinetics, there is no a priori reason why the barrier 
entropy cannot become negative, however the discussion 
in sec. IIIBI in terms of extreme value statistics, indi- 
cates the analogy to thermodynamics may be also applied 
to the phenomenon of freezing. This has been demon- 
strated by Bouchard and Mezard— where the freezing 
phenomenon could be understood from the perspective 
of extreme value statistics. These aspects are also briefly 
considered in appendix C. 

Alternatively, the entire procedure entailed by eqns. lA3l 
to IA6I can be viewed as a change of integration variables 
resulting in the integral of eqn. |9]in sec. IIIBI being trans- 
formed to a contour integral in the complex plane. This 
approach becomes advantageous when considering the 
modified Gaussian of sec. Ill Dl which gives the barrier 
entropy as 



S{E) = kBN 



_ l(g{E)/N- 



SE 




(All) 



and the relation between internal barrier energy and tem- 
perature as 



g{E)/N~E 

se' 



9\E) 



E=U{T) 



f (A12) 



To obtain the internal energy as a function of temper- 
ature, and thereby be able to construct the free barrier 
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energy as a function of temperature, eqn. IA12I must be 
inverted. For particular forms of g{E) this may be done 
analytically, but more generally, such an inversion is done 
numerically. For the case of the barrier energy at freez- 
ing, S{Ef) — leads to 

Ef = Ng-^ (e - SEV2aj . = Ng-^ (sf^j (A13) 

enei 

pure Gaussian distribution (eqn. [2T|). 



t( 



Q 

where Ef is the freezing barrier energy derived from a 



Appendix B: A simple deformation model 

To gain an estimate of the relationship between Toxp 
and the strain rate £totai the following simplified vis- 
coplastic model is developed. Assuming an additive elas- 
tic and plastic strain rate, the total strain rate is given 
by 



Etotal = Ec +ep(o-) = ^ +£p(o-), (Bl) 

which gives the first order equation 
da 



— = G [etotai - £p(o')] 



(B2) 



whose solution may be written as an integral equation 

ct+St 



a{t + St) = cr{t) + G 



dtip{a{t)) 



(B3) 



For a small enough time interval this may be approxi- 
mated as 



a{t + St) ^ a{t) + Gitote^iSt 



1 



£total 



(B4) 



which can be iterated to give a, constant strain-rate, 
stress-strain curve. Inspection of the above equation 
gives the condition for perfect plastic flow {a{t + St) = 
(j{t)) as when the plastic strain rate is equal to the total 
applied strain rate (£p = etotai)- 

To obtain an estimate for ip, and therefore a connec- 
tion to £totai, the simplest approximation is to linearly 
relate the plastic strain rate of a glass of volume V to 
the presently derived plastic rate, [Tp]~^: 



Ep = Ae X [TpY 



V 



N{lo) 



(B5) 



In the above, Iq is a measure of the average atomic spac- 
ing, giving V/N{lo)^ as the number of heterogeneous vol- 
umes within the material volume V, and Ae is an esti- 
mate of the characteristic plastic strain arising from an 
irreversible structural transformation occurring within a 
heterogeneous volume. Under the assumption that an 
irreversible structural transformation can be well repre- 
sented by a localized shearing of the material, according 



Ae will be given by the ratio of the characteristic 
slipped area projected onto the gauge cross-section mul- 
tiplied by the slip distance, and V. Assuming the charac- 
teristic slipped area has a length scale that is comparable 
to the heterogeneity length scale (say, half) and the cor- 
responding slip distance will be of the order of an atomic 
spacing (since the mediating /3-relaxation processes will 
involve 0(1) atomsSI), the above equation reduces to 



1 



^total 



4]Vl/3 



(B6) 



For Vitreloy-1 the proportionality constant is approxi- 
mately 0.01 when TV ~ 3000 (corresponding to a « 0.1). 
It is noted that all geometric factors have been ignored 
and that the above should be viewed at best as an order 
of magnitude estimate, which is sufficient given that the 
resulting timescale enters enters the logarithmic factor 
in eqn. 1361 In the present context, a choice of a corre- 
sponding characteristic plastic strain, Ae, is actually a 
choice on the nature of the mega-basin the material fi- 
nally enters upon exiting its current mega-basin — an 
aspect that is not included in the current model. The 
current choice simply assumes that the system finds (on 
average) a final state that is compatible with the external 
loading geometry. 



Appendix C: The connection to extreme value 
statistics 



When sampling the barrier energy distribution P{E), 
M times, the probability that the minimum barrier en- 
ergy is i?* is given by2^ 



PnuniE*) = MP{E*) {1-P<{E*)) 



M-1 



(CI) 



where P< (E) is the repartition probability or the cumu- 
lative distribution function of P{E): 



P<{E*) 



E* 



dEP{E). 



(C2) 



Eqn. lClj is the required extreme value distribution associ- 
ated with the statistics of the smallest barrier energy in a 
heterogeneous volume containing M possible barrier en- 
ergies. Fig. [3]plots eqn. ICll using P{E) derived from the 
parameters of Vitreloy-1 at a number of different stresses 
spanning both the low and high temperature plastic tran- 
sition stress versus temperature curve of fig.|51 Inspection 
of this figure reveals distributions that are peaked at en- 
ergies close to that of the kinetic freezing barrier energy 
scale and with increasing stress are shifted to smaller bar- 
rier energies. The shape of the extreme value distribution 
also changes upon increasing the stress, where for the 
largest stress the distribution could be fitted optimally 
using a Weibull distribution, and for the lowest stress, a 
Gumbel distribution provided the best fit. Such a result 
is expected since at low stresses (the high temperature 
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large M, 



E (eV) 



FIG. 9. Plot of extreme value barrier energy distribution de- 
rived from the barrier energy distribution for the parameter 
set corresponding to Vitreloy-1 (fig-Et) for five stresses span- 
ning both the low and high temperature regimes (fig. [8|). 
Also included are the optimal Gumbel fit for the lowest stress 
distribution and the optimal WeibuU fit for the highest stress 
distribution. 



regime) the Gaussian form dominates with the extreme 
value statistics unaware of the necessity of a positive bar- 
rier energy, giving a Gumbel extreme value distribution. 
However as the stress increases (and the low temperature 
regime is entered) , the modified Gaussian form begins to 
dominate resulting in a WeibuU extreme value distribu- 
tion. 

Inspection of the corresponding mean derived from 
each distribution reveals values that are consistently 
higher than that predicted by the kinetic freezing bar- 
rier energy eqn IA13I The origin of this may be seen by 
determining an approximate expression for the extreme 
value distribution derived mean value. For sufficiently 



(l-P<(i?*))*'^^exp(-AfP<(i?*)) 



(C3) 



The distribution entailed by eqn. ICll will peak when the 
argument of the exponetial in the above is approximately 
one giving a definition for the mean minimum value, , 
as 



p<(i?;) 



1 



(C4) 



For small enough E, the repartition distribution can be 
approximated to logarithmic accuracy as 



P<iE) 



EP{E) 



27r 



(C5) 



resulting in eqn. IC4I reducing to 



iog(p(ii;;))~iog 



(C6) 



For the pure Gaussian form, eqn. [51 this reduces to 




(C7) 



The above equation gives an accurate estimation of the 
average minimum energy barrier, E^ , for the extreme 
value distributions for Vitreloy-1 shown in fig. [HI It is 
only in the bulk limit, N oo, that eqn. IC7I becomes 
identical to eqn. IA8I resulting in Ef = Ef and a for- 
mal equivalence between extreme value statistics and the 
thermodynamic treatment of the kinetic freezing phe- 
nomenon developed in appendix A, as was shown by 
Bouchaud and Mezard^. The right-hand-side of eqn. IC7I 
is therefore a correction that takes into account the finite 
size of the heterogeneous volume. 
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